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The  bulk  of  the  research'  performed  under  this  contract  has 
dealt  with  mathematical  problems  of  numerical  ocean  acoustics,  of 
tiie  kind  described  in  Ref.  1,  pp.  64  3- 54 .  These  concerned  the 
propagation  of  sound  waves  in  (generally  inhomogeneous)  elastic 
fluids,  with  special  reference  to  the  consistency  of  the  elastic- 
fluid  model  with  'ray  theory'  (Fermat -Huygens) ,  in  predicting 
ref  1  ect  i  on  ,  ref  raft  i  on  ,  and  d  i  f  f  ract.  ion  ■  vs  '  ('**,' 

I  have  written  up  most  of  the  fruits  of  this  research  for  pub¬ 
lication  in  Refs.  2-4,  of  which  Ref.  2  lias  appeared  (reprint  en¬ 
closed).  A  copy  of  Ref.  3  has  been  retyped  at  t lie  Naval  Underwater 
Systems  Center  for  photographic  reproduction.  Typescripts  of  Refs. 
4  and  5  have  been  sent  to  the  editors  of  the  Proceed  lugs  in  which 
they  will  appear.  Much  of  Ref.  4  is  a  Jess  detailed  but  more 
polished  exposition  of  results  presented  in  Refs.  2  and  3.  These 
results  include  the  following: 

A.  Formulation  of  a  precise  mathematical  definition  of  the  concept 
of  an  inhomogeneous  el  a  s  t i c  f 1 u i d  (Rei .  3,  Sec.  9  and  Appendix 
C ;  cl.  Ref.  2,  Secs.  10  and  IS)  . 

B.  Derivation  plane  of  exact  second-order  linear  d i \  ! e ren t  I  a  1  cqo.u 
lions  for  waves  in  i nhomogeneous  Chaplygin  fluids.  (Rei.  3, 
Secs.  JR  end  14;  Ref.  !,  Sec.  0). 

C.  The  conclusion  that,  in  more  than  one  d  i  men  s  i  on  .  sound  waves  '  u 
such  fluids  will  usually  generate  first-order  vert  icily  (Ref.  3 
Secs.  9  and  11).  Thus  the  umk;  1  a sumpt  i  on  o:  the  oxi.len'e  e  : 
a  velocity  potential  is  not  justified,  c\en  if  there  is  r.o 

v  i  s  c  o  s  i  t  y  . 

I)  'the  tent  at  ivc  conclusion  that  ,  novr  i  t  be  I  ess  ,  i  os  l  '  and 
Huv.’eus  '  l’i  i:u  i  p  1  e  s  a  re  asympl  oi  it;'  !  \ >  a  i  id  ‘  .  \  o  ry  sho-'i 
wave:-  in  a  ri>  i  l  ra  l  y  ii  astir  111!  ids  t  !'.  ■  i' .  1,  u  )  . 
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E.  Clarification  of  the  physical  significance  of  the  Sommcrfeld- 
Runge  cikonal  function  for  sound  waves  in  fluids  (Ref.  2,  Sec. 
10;  Ref.  3,  Sec.  5;  Ref.  4,  Sec.  6). 

Namely,  each  "eikonal  surface"  on  which  the  Sommer f el d ~ Runge  eikonal 
function  assumes  a  constant  value  is  asymptotically  an  i sohar i c 
surface  of  constant  pressure  at  all  times,  associated  with  "simple" 
sound  waves  (Ref.  1,  p.  64  6]  in  which  the  fluid  is  constrained  work- 
lessly  to  move  along  acoustic  rays. 

The  preceding  results  hold  in  static  fluids.  In  mo v i ng  fluids 
(Ref.  3,  Secs.  7-8),  although  ray  theory  can  still  be  applied,  the 
use  of  variational  principles  involves  tine  concept  of  a  F i n s  1  e r 
space ,  and  very  unfamiliar  geometrical  ideas. 

The  elastic  fluid  models  discussed  above  all  conserve  mechani¬ 
cal  energy;  they  do  not  predict  any  absorption  or  d  i  spurs  i.  on .  I 
have  made  a  critical  analysis,  in  historical  perspective,  of  several 

standard  treatises  concerned  with  models  intended  to  explain  there 

*  * 

related  phenomena.  My  impressions,  reported  in  Ref.  3  (Secs.  3-S 
and  Appendix  A]  and  Ref.  4  (Secs.  7-11),  arc  summarized  in  the 
Appendix  to  the  present  Report.  1  have  correlated  them  with  the 
discussion  of  absorption  and  dispersion  (primarily  in  gases)  on  pa. 
5>49-5>62  of  Allan  Pierce's  admirable  back  Acoustics. 

My  main  conclusion  is  that  the  sinnda.rd  modern  explanation  in 
terms  of  re  lax. at  ijaii^t  ines,  although  six')-  years  old,  has  not  tel  In  e 
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.substantiated  (especially  in  liquids)  by  clear  answers  to  many  b;si 
questions.  These  include  the  following: 

F.  To  what  extent  is  the  absorption  of  sound  per  wave  length,  c-.,\  , 

in  air,  CO),  and  other  dilute  gases  determined  by  the  absolute 
temperature,  T,  and  the  ration  f/p  of  the  frequency  to  the 
pressure?  Vj 

G.  To  what  extent  are  contributions  to  a  from  different  causes 
demonstrably  additive ,  in  gases  and  in  liquids? 

H.  How  well  can  one  predict  the  locations,  breadths,  and  heights  o 
the  two  graphs  drawn  in  Fig.  10-12  of  Pierce?  How  well  are 
these  correlated  experimentally,  in  gases  and  in  liquids? 

k,  How  are  the  hu  1  k  viscosities  of  liquids  and  gases  best  defined 
and  measured?.  The  old  book  by  A.  Harr,  A  Manual  of  Viscometrv, 
on  shear  viscosity,  sets  a  good  example. 

Compu t a t i on a  1  a c ou s tics.  When  this  contract  was  set  up,  I 
planned  to  polish  my  1983  soft- cover  monograph  N’umeri  ca  1  )■' lu:i  d 
Dynami  cs  for  publication  in  hard-cover  form,  including  in  it.  new 
material.  A  paper  submitted  to  SIAM  Review  on  "Difference  methods 
for  solving  convect ion-dif fusion  equations,"  ccauthored  by  E.C. 
Garland,  or.,  and  R.F.  Lynch,  was  written  with  this  in  mind.  So  wa 
t  tie  visit  to  Harvard  by  Prof.  Lynch  in  the  summer  of  1 9  S  6 ,  during 
wtii  eh  we  discussed  improving  difference  approx  i  mat  ions  to  Neuman; 
and  mixed  boundary  conditions.  Since  then.  Prof.  Lynch  ha^  inciwp- 
oraied  those  into  the  Purdue  Fit. I, PACK  package ,  and  written  two  paper 
on  t  ho  sub j  ee l . 

1  planned  La  do  extensive  rewriting  in  om'-nt  ion  with  a  Svi„in 
on  ,\mi!i  i  i  r.  I  i-luid  Dynamics  which  i  was  schedi.Med  1  a  lead  at  t  he 
N.\al  l\m  i  yradaa  !  e  Sc  hoe  I  laM  winter.  Ilo.vev.  i.  ini'-  piojt.cl  p :  -  . 
res.-.ed  s  I  :*w  >  y ,  and  aider  a  heart  at  t  .m  L  made  in  -  •>  i  ve  up  ibis  So 
:  ii  l!ec(  ;  .lev  ,  I  abandoned  it.  Instead,  I  sju  in  li  e  pm  imrin:;  Re 1  . 
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over  half  of  which  is  devoted  to  publicizing  the  major  role  played 
by  navy- sponsored  research  in  the  crucial  early  development  of 
scientific  computing.  This  invited  paper  will  appear  in  a  volume 
which,  I  believe,  will  be  widely  read. 


Appendix.  Absorption  and  Dispersion 

Already  by  1860,  general  mathematical  theories  of  viscosity 
(Stokes)  and  heat  conduction  (Fourier)  had  been  constructed.  Both 
of  these  theories  predict  that  sound  waves  in  fluids  should  be  sub¬ 
jected  to  absorption  and  dispersion,  and  that  the  rate  of  internal 

2 

sound  absorption,  «,  should  be  proportional  to  f  ,  the  square1  of 
the  frequency.  Actually,  the  two  theories  can  be  incorporated  into 
a  more  elaborate  theory  of  homogeneous  vi scorns  thermo el ast  ic  fluids. 
This  theory  involves  three  parameters  besides  those  which  character¬ 
ize  elastic  fluids:  the  thermal  conductivity  k,  the  shear  viscosity 
p,  and  the  bulk  viscosity  p ’ . 

Until  1925,  this  theory  seemed  adequate,  partly  because  the 
internal  absorption  and  dispersion  of  sound  in  air  and  water,  on  a 
laboratory  scale  (<  10m),  are  almost  imperceptible  over  the  audible 
range  (roughly. 25  Iiz-5  kllr.)  .  The  effects  of  t.  and  n'  were  simply 
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waves  in  the  range  5  kHz-]  milz  revealed  rates  of  internal  sound 
absorption  in  fluids  of  Ilia  -  100a^  or  even  more.  These  are 
generally  attributed  to  molecular  relaxation  effects,  associated 
with  "relaxation  times"  t  required  to  transfer  "rotational"  to 
"translational"  energy.  In  gases ,  this  concept  can  be  rationalized 
in  terms  of  Maxwell -Boltzmann  mass-spring  model  s  of  molecules;  more 
over,  kinetic  theory  suggests  that  Aa  should  be  a  single-valued 
function  of  the  ratio  f/p  of  the  sound  wave  frequency  to  the  gas 
pressure.  From  1925  to  1940,  discoveries  of  major  deviations  from 
(1)  in  gases  and  liquids  were  correlated  with  developments  in 
theoretical  molecular  physics,  which  was  being  revolutionized  by 
quantum  mechanics. 

After  World  War  II,  intense  interest  in  shock  waves  helped  to 
stimulate  extensive  further  studies  of  the  absorpt i on  of  mechanics 
energy  in  80  or  more  chemically  pure  fluids.  Liquids  were  class  i 
fied  into  qualitatively  different  types,  as  regards  their  sound 
absorption  properties.  In  addition,  important  effects  of  inipuj  iti 
and  sound  absorption  i  a  mixture's  (such  as  air)  continued  to  be 


st ud i ed  . 

Bulk  \iscosity  also  received  careful  scrutiny  for  the  first 
time  after  1  940.  'Ihus  i, .  1  .  Taylor  (Scicntilic  Papers,  iv,  '"..’3) 
showed  that  air  bubbles  in  water  can,  theoretically,  make  !•'  :  <» V <»*•'.  . 

(Hi  fiord  Trues do  1  1  made  the  lirst  rigorous  ana  1  y  s  i  :•  ol  the  theeie.  , 
e  a  1  depen.lei.ee  of  u  on  i  I  e  d  i  mens  i  on  1  'ess  rat  i  os  end  l'r  !  |.  ■ 

in  a.  gene-rs  I  \  i  some  t  he  nw  !  a  s  t  ic  fluid,  and  of  '  hoo  ret  i  a  1  dev  i .  i  - 
1  )  .  1!  i  s  c  i  i  t  i  quo  of  p  r.i.  ed "  r>  mod  In  lit  " 


lions  f i or  l 


need  1 o  fit 


l  e  I  s  :■  e  i  u 


times"  to  empirical  data  has  never  been  clearly  answered.  It  would 
be  desirable  to  have  a  clearer  idea  of  the  magnitude  of  Trucsdcll's 
corcctions  to  formula  (1),  which  is  still  believed  to  hold  (with 
v  '  =  0)  for  inert  gases  like  He. 

Of  even  greater  interest  is  the  empirical  evidence  for  Pierce’ 
f ormula  (10-8.11) , 

(2)  a(f)  =  ac£(f)  +  l  av  (f)  , 

where  different  dominate  over  different  ranges  of  f,  and  his 

analogous  formula  (10-8.16)  for  dispersion.  It  is  not  clear  for 
which  gases  and  liquids  these  formulas  have  been  proven  experiment¬ 
ally  to  hold,  with  what  precision  over  what  ranges  of  f  and  state 
(p,T).  Likewise,  for  which  pure  and  impure  fluids  have  the  param¬ 
eters  (a  A)  and  t  which  occur  in  these  formulas  been  deduced 
v  in  v 

from  general  basic  formulas  of  quantum  mechanics,  and  with  what 
precision? 

Until  such  questions  (which  are  related  to  Questions  P-l 
stated  in  the  body  of  this  report)  are  answered,  1  think  that  tin.' 
status  of  ’relaxation  time’  rat iona 1 i cut  ions  of  absorption  and  dis¬ 
persion  phenomena  resembles  that  of  ’classical’  rational  i  /.at  ions  as 
of  1  POD .  Although  highly  plausible,  their-  range  of  validity 
(especially  for  liquids)  is  unclear! 
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